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An infectious bronchitis virus (IBV), ck/CH/LZJ/111113, was isolated from a HI 20- 
vaccinated chicken which showed disease suspected of IBV infection. Neutralization 
testing showed that ck/CH/LZJ/111113 was distinct from either the Chinese predominant 
IBV LX4-type or Mass-type vaccine strains. Phylogenetic analysis confirmed that ck/CH/ 
LZJ/111113 is of the 4/91 type; however, further extensive analyses of full-length genomes 
identified occurrence of recombination events. Therefore, ck/CH/LZJ/111113 originated 
from the recombination events between ck/CH/LDL/091022- and 4/91-like strains at three 
switch sites located upstream of the spike (S) glycoprotein gene, and the 3' ends of SI and 
nuceocapsid (N) genes, respectively. The difference of serotypes and tissue tropisms in 
kidneys between ck/CH/LZJ/111113 and ck/CH/LDL/091022 may have been contributed by 
the uptake of the SI gene by a ck/CH/LDL/091022-like virus from a 4/91-like strain. This 
recombination event took place at the 3 ' end of the N gene and the 3 ' untranslated region 
may account for differences in replication efficiency in tissues of chickens inoculated by 
the two viruses. 

© 2012 Elsevier B.V. All rights reserved. 


1. Introduction 

Coronaviruses (CoVs) are a family of enveloped 
positive-stranded RNA viruses and found in a wide variety 
of animals in which they can cause respiratory, enteric, 
hepatic and neurological diseases of varying severities. 
Based on genotypical and serological characterizations, 
CoVs were classified into three genera, Alphacoronacirus, 
Betacoronavirus and Gammacoronavirus (Carstens, 2010). 
Avian infectious bronchitis virus (IBV) belongs to Gamma- 
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coronavirus genus and causes a highly contagious upper 
respiratory disease, a decrease in egg production, poor 
eggshell quality, reduced hatchability, nephritis and 
sometimes false layers in chickens. 

CoVs contain an unsegmented, single-stranded, posi¬ 
tive-sense RNA genome of 28-32 kb, with a 5' cap and a 3' 
polyA tail. In the IBV genome, similar to other CoVs, two- 
thirds of the 5' region encodes the la and lab polyproteins, 
which are proteolytically cleaved by two virus-encoded 
replicase proteins (the papain-like and 3C-like protei- 
nases) into at least 15 nonstructural proteins (nsp2- 
nspl 6) (Ziebuhr et al., 2001 ). In the remaining one-third of 
the genome, most coronaviruses, including IBV, code four 
major structural proteins; the spike (S) glycoprotein, the 
small envelope (E) protein, the membrane (M) glycopro¬ 
tein, and the nucleocapsid (N) protein. It is generally 
believed that the SI subunit contains the receptor-binding 
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activity of the S protein, carries virus-neutralizing and 
serotype-specific determinants, and is a determinant of 
cell tropism. Interspersed among the structural protein 
genes are two small accessory protein genes (genes 3 and 
5) that vary in number and sequences among the IBVs. The 
5' and 3' untranslated regions (UTRs) usually harbor 
important structural elements, which are involved in 
replication and/or translation (Masters, 2006). 

LX4 (QX-like) and 4/91 are two of the most well 
characterized and important IBV serotypes that were first 
isolated in 1998 in China (Liu and Kong, 2004) and 1985 in 
France (Cook et al., 2012), respectively. Viruses of these 
serotypes have been detected in several European as well 
as Asian countries, including China, although at different 
frequencies (Liu and Kong, 2004; Cook et al., 2012), causing 
heavy economic losses. Here, we isolated and character¬ 
ized an IBV strain (ck/CH/LZJ/111113) and reported data 
from antigenic and genomic analyses and replication 
efficiency in chicken tissues and showed that it represents 
a recombinant coronavirus infectious bronchitis virus. 

2. Materials and methods 

2.2. Viral isolation 

In the course of our continuous surveillance activities 
for IBV in China in 2011, an IBV strain, designed as ck/CH/ 
LZJ/111113, was isolated from a diseased trachea sus¬ 
pected of an IBV infection. The chicken flocks were 
vaccinated against IB with commercial live attenuated 
HI 20 vaccine and showed early signs of respiratory disease 
at 20 days old. The gross examination showed mild to 
severe tracheitis and proventriculitis. 

Viral isolation was performed as previously described 
(Liu and Kong, 2004). Blind passages were performed until 
the characteristic embryo changes appeared, such as 
dwarfing, stunting, or curling of the embryos, which were 
observed between 2 and 7 days post-inoculation, based on 
a previous report. In addition, the IBV strain ck/CH/LDL/ 
091022 (Sun et al., 2011) was used in this study. The 
embryo-propagated viral stocks of ck/CH/LDL/091022 and 
ck/CH/LZJ/111113 were produced by inoculating the virus 
into embryonated SPF chicken eggs via the allantoic cavity 
and collecting the infectious allantoic fluid 72 h post¬ 
inoculation as previously described (Liu and Kong, 2004). 

2.2. Eggs and chicks 

Fertile White Leghorn SPF chicken eggs were obtained 
from the Harbin Veterinary Research Institute, as were the 
White Leghorn SPF chicks. The birds were maintained in 
isolators with negative pressure, and food and water were 
provided ad libitum. 

2.3. Virus-neutralization tests 

The titers of the two viruses were determined by 
inoculation at 10-fold dilutions into groups of five 10-day- 
old embryonated chicken eggs. The 50% embryo infectious 
dose (EID 50 ) was calculated by the method of Reed and 
Muench (1938). The (3 VN method with constant virus and 


diluted serum was employed in SPF chicken embryos for 
serotyping. The R-values were calculated by the formula of 
Archetti and Horsfall (1950). 

2.4. Sequencing and phylogenic analysis of IBV strains ck/CH/ 
LDL/091022 and ck/CH/LZJ/111113 

For complete genomic sequencing, viral RNA was 
extracted from 400 |ul1 of ck/CH/LDL/091022 and ck/CH/ 
LZJ/111113 viral stocks, respectively, as described pre¬ 
viously (Liu et al., in press), using TRIzol reagent according 
to the manufacturer’s instructions. Fourteen overlapped 
PCR fragments spanning the entire viral genome were 
amplified using specific primer sets (Liu et al., in press). The 
3' and 5' ends of the viral genomes were confirmed by 
rapid amplification of cDNA ends using a 3'/5' RACE kit 
(TaKaRa) according to the manufacturer’s instructions. The 
sense and antisense primers used for amplifying the 3' and 
5' ends of the genome had been designed as previously 
described (Liu et al., in press). The PCR products were 
cloned into the pMD-18T vector (TaKaRa) following the 
manufacturer’s instructions, and each fragment of the viral 
genome was sequenced at least five times and the 
consensus sequence was determined. Prediction of the 
ORFs was conducted using Vector NTI Advanced 10 
(Invitrogen), and the sequences were analyzed using 
Lasergene DNAStar (version 7, Lasergene Corp, Madison, 
WI, USA). The pairwise nucleotide identity was determined 
using Vector NTI Advanced 10 and multiple sequence 
alignments were generated using Clustal-W (Thompson 
et al., 1997). 

A total of 30 or 39 IBV reference strains for which the SI 
or entire genomic sequences were available in the 
GenBank database were selected for phylogenetic analysis. 
The selected avian coronavirus reference strains and their 
accession numbers are shown in Fig. 1A and B, respec¬ 
tively. 

2.5. Accession number 

The genomic sequences of the IBV strains ck/CH/LDL/ 
091022 and ck/CH/LZJ/111113 have been submitted to the 
GenBank database and have been assigned the accession 
numbers JX195175 andJX195176, respectively. 

2.6. Experimental design 

Ninety, 1-day-old SPF White Leghorn chicks were 
housed in different isolators and divided into three 
groups, each containing 30 birds. Chickens in groups 1 
and 2 were inoculated with either ck/CH/LDL/091022 or 
ck/CH/LZJ/111113 (lx 10 s EID 50 per chick, respec- 
tively) by oculonasal administration at 15 days of age 
with a dose of loglO 6 EID 50 per chick. Birds in group 3 
were mock-inoculated with sterile allantoic fluid and 
served as the control. Ten birds from each group were 
killed humanely at 4 and 7 days post-inoculation. 
Tissue samples were collected from the trachea, lung, 
kidney, Harderian gland, thymus, spleen, bursa of 
Fabricius, cecal tonsil, proventriculus, duodenum, liver, 
pancreas, small intestine, large intestine, rectum and 
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Fig. 1. Relationships among the ck/CH/LZJ/111113, ck/CH/LDL/091022 and reference strains illustrated by a maximum likelihood phylogeny unrooted tree, 
based on SI nucleotides 1-1626 (A) and the full-length of genomic sequences (B), respectively.In (A), the strains ck/CH/LZJ/111113, ck/CH/LDL/091022 are 
shown in bold. The genotypes of 4/91-, LX4- and mass-type strains were shown in pink, blue and red, respectively. In (B), the US and China strains were 
shown in pink and blue, respectively. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of the 
article.) 


caecum for viral RNA detection by real-time RT-PCR. 
The remaining birds in each group were examined daily 
for signs of infection for 30 days post-inoculation. Blood 
samples were collected from all birds in each treatment 
group at 4, 8, 12, 16, 20 and 24 days post-inoculation. 
Serum samples were assayed using a commercial total 
antibody ELISA (IDEXX Corporation, Westbrook, ME, 
USA) according to the manufacturer’s instructions. Each 
sample was tested in triplicate. 

2.7. Real-time RT-PCR 

The homogenized tissue samples were used for 
quantitative real-time RT-PCR as described previously 
(Jones et al., 2011). One-step real-time RT-PCR reactions 
were performed using One Step PrimeScript® RT-PCR kit 
(TaKaRa) on a LightCycler® 480 real-time PCR system 
(Roche, Basel, Switzerland). The data were analyzed using 
LightCycler® 480 Software Version 1.5. All real-time RT- 
PCR products were confirmed by electrophoresis on a 2.0% 
agarose gel and visualized using ethidium bromide. 

3. Results 

3.2. Serotype of IBV by cross virus neutralization tests 

The ck/CH/LZJ/111113 was serologically characterized 
with the Chinese predominant LX4-type strain, ck/CH/LDL/ 
091022, and the vaccine strain, HI20. Cross virus 
neutralization tests were performed using anti-sera 


against strains ck/CH/LZJ/111113, ck/CH/LDL/091022 and 
HI20 to determine the antigenic relationship among the 
three strains. End-points were calculated by Reed and 
Muench method (1938). The R value for strains ck/CH/LZJ/ 
111113 and ck/CH/LDL/091022, and ck/CH/LZJ/111113 
and HI20 were 28.42 and 29.73%, respectively, which 
indicated that ck/CH/LZJ/111113 is distantly related to 
strains ck/CH/LDL/091022 and HI20. Thus, strain ck/CH/ 
LZJ/111113 represents a new serotype that is antigenically 
distinct from LX4 and Mass serotypes of IBVs. 

3.2. Phylogenetic analysis of the strain ck/CH/LZJ/l 11113 

The SI gene of ck/CH/LZJ/111113 had 99.6 and 99.5% 
identity with pathogenic and attenuated UK/4/91 strains, 
respectively (Lig. 1 A). These figures are very similar to the 
comparison of a recent Indian 4/91-type isolate, India/ 
NMK/72/IVRI/10. Thus, the SI sequence of ck/CH/LZJ/ 
111113 was clearly of the 4/91 type. However, the 
sequence was markedly different from those of 4/91 from 
position 1453 in the SI gene to the end of the SI gene, from 
which position it was closely related to a LX4-type virus, 
ck/CH/LDL/091022. 

The ck/CH/LZJ/111113 and ck/CH/LDL/091022 
sequences were assembled into one contiguous sequence 
to represent the entire viral genomes. The sequences of 
27,675 and 27,673 nucleotides were obtained from ck/CH/ 
LZJ/111113 and ck/CH/LDL/091022, respectively, including 
the poly-A tail at the 3' end. The genomes of both viruses 
were similar overall in their coding capacities and genomic 
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A. I 


4/91 ( FN811147) 

ck/CH/LZJ/111113 
ck/CH/LDL/091022 


TAT TTCGG 

GCAGCGTCTTTGGTGCTTGCCCGTAAACATACTAATTGTTGCACATGGTCTGAGCGCATTTATAGATTGT 14303 

14342 


4/91 (FN811147) 

ck/CH/LZJ/111113 
ck/CH/LDL/091022 


T TGTAG GT G 

ACAATGAATGCGCCCAAGTATTGTCTGAAA...AACCAGACTTAAAAAAGGGCCCACATGAATTTTGTTC 

C T 


4/91 (FN811147) 

ck/CH/LZJ/111113 
ck/CH/LDL/091022 


ACAGT TAAG T C 

ACAGCACACAATGCTTGTAGAGGTGGATGGTGAGCCTAGGTACTTACCATACCCAGACCCCTCACGCATT 

A T C A A G A 


14848 

14887 


A. II 

4/91(JN192154) - NA - 

ck/CH/LZJ/11113 TTGTTAATTTGAAAACTGAACAAAAGACAGACTTAGTCTTTAATTTAATTAAGTGTGGTAAGTTACTGGTAAGAG 20367 
ck/CH/LDL/09102 A C TG C C A G AAC C TAT 20416 


*** (Initiation codon of spike gene) 

4/91(JN19 215 4) 

ck/CH/LZJ/11113 ATGTTGGGCAAACCGCTTTTACTAGTGACTCTTTGGTATGCACTATGTAGTGC-TTTGCTTTATGATAAAAAT 

ck/CH/LDL/09102 G GT A G TT CA T TG AAA T G TC GCC 


4/91(JN19 215 4) 
ck/CH/LZJ/11113 
ck/CH/LDL/09102 


T 

ACTTACGTTTACTACTACCAAAGTGCCTTTAGG...TTTGTTGTTCGAGGTGCA TATGGTCTTAATTATTA CAAGG 
A T G A AG CAGC 


4/91(JN192154) 
ck/CH/LZJ/11113 
ck/CH/LDL/09102 


T A C C T C A CAT 

TCAGTCCTTGTGAAGATGTTAACCAACAGTTTGTAGTGTCTGGTGGTAATATAGTTGGTATTCTTACTTCTAGAAA 
A C C 


4/91(JN192154) A TTTA CA C AC TCA T AC AA 

ck/CH/LZJ/11113 TGAAACAGGTTCTGAACAGGTTGAGAACCAGTTTTATGTTAAGTTAACCAATAGCTCACATCGTCGCAGGCGTTCT 
ck/CH/LDL/09102 

(cleavage recognition site)############## 

4/91(JN192154) G ACTGGG T A CC C T G 

ck/CH/LZJ/11113 ATTGGCCAAAATGTAACAAGTTGTTCTTATGTT...GTTGTGGATGCTTTGGTATCATTCCTCTAATGAGTAAGTG 
ck/CH/LDL/09102 


4/91(JN192154) T C A 

ck/CH/LZJ/11113 TGGTAAGAAATCATCTTATTACACGACTTTTGATAATGATGTGGTGACTGAACAATACAGACCTAAAAAGTCTGTT 23860 
ck/CH/LDL/09102 23911 


A. Ill 


*** (Initiation codon) 

4/91 (GQ227602) -GA A CG A AGAATAT 

ck/CH/LZJ/111113 ATGTCCAACGAGACAAATTGCACCCTTGATTCTGAGCAGGCTATTCTGCTCTTCAAGGAATATAATCTATTTATAAC 24585 
ck/CH/LDL/091022 24636 


4/91 (GQ227602) T GT C C 

ck/CH/LZJ/111113 CGCATTCCTACTGTTTATAACCATACTACTTC...ATACTGGTGACCAAAGCGGAAATAAGAAAAGGTTTGCTACAT 
ck/CH/LDL/091022 


4/91 (GQ227602) 
ck/CH/LZJ/111113 
ck/CH/LDL/091022 


G T T 

TTGTCTATGCAAAGCAGTCAGTAGACACTGGCGAGCTAGAAAGTGTAGCAACAGGAGGAAGTAGTCTTTACACATAA 25186 

(Stop codon)*** 25237 


Fig. 2. Analysis of recombination among IBV ck/CH/LDL/091022, 4/91 and ck/CH/LZJ/111113 strains. (A) Multiple sequence alignments surrounding 
numbered regions in A. The numbers on the right of each alignment showed the nucleotide positions in genome of each virus. The sequences of ck/CH/LZJ/ 
111113 are listed and the only nucleotides differing from those of ck/CH/LZJ/111113 are depicted. The region where the template switches (switch site) 
have taken place is underlined. The deleted nucleotides are indicated by -. NA represents unavailable sequence. (B) Percentages of nucleotide sequence 
identity among ck/CH/LZJ/111113, ck/CH/LDL/091022 and 4/91.Percentages of nucleotide sequence identity of gene fragments before and after the switch 
sites in SI and N genes are also indicated. 
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A. IV 


4/91(EU780081) 
ck/CH/LZJ/111113 
ck/CH/LDL/091022 


***(Initiation codon of N gene) 

A GCG AAGA G CAA 

ATGGCGAGCGGTAAAGTATCTGGAAAGTCAGACTCCCCCGCGCCAATCATCAAACTAGGAGGGCCTAAAC 26004 

A 26055 


4/91(EU780081) 
ck/CH/LZJ/111113 
ck/CH/LDL/091022 


A 


A 


A 


G A 


C 


G C 


CACCAAAGGTAGGGTCATCTGGAAGTGCA...GCCCAAAAGAGTAAAGATTTGTGATGAGTGTATTGATG 


A 


G 


4/91 (EU780081) A G T T AACCG GCA A 

ck/CH/LZJ/111113 GTGTAGGCACACGCCCAAAAGACGAAGTTGTAAGACCAAAATCACGCTCAAGTTCAAGACCTGCT ACAAG 
ck/CH/LDL/091022 A T 


4/91(EU780081) 

ck/CH/LZJ/111113 AGGAAATTCTCC TGCGCCAAGACAACAGCGCCAAAAGAAGGAGAAAAAGCCAAAGAAGCAGGAAGATGAT 
ck/CH/LDL/091022 A A G A 

4/91(EU780081) 

ck/CH/LZJ/111113 GTAGATAAGGCATTGACCTCAGATGAAGAGAGAAGGCATTGACCTCAGATGAAGAGAGGAACAATGCACA 
ck/CH/LDL/091022 G A G A G 

4/91(EU780081) 

ck/CH/LZJ/111113 GCTCGAGTTTGATGATGAACCCAAAGTGATTAATTGGGGAGATTCAGCTCTTGGGGAAAATGAATTGTAA 27164 
ck/CH/LDL/091022 AA TAT CTG 27216 

(Stop codon of N gene)*** 


A. V 

4/91(HM210804) 

ck/CH/LZJ/111113 GTAACATAATGGACTTGCTGCATTTGCTGTCACATTTTGTTAAATACTATTTTGTTGCTTTCCTATCAAT 27234 
ck/CH/LDL/091022 A TA GG GCTG TGT T C CTC T C 27285 

4/91(HM210804) 

ck/CH/LZJ/111113 TATTACAGGCATTGATTGTGATTATGTGCAAT...TATTGTGATTCTCATTAGTTTGCTTTATCGTAGAA 
ck/CH/LDL/091022 C TG TT G G T A TC G TTA TC T A 


4/91(HM210804) 

ck/CH/LZ J/111113 GTGCAATAGTAAGAGTTAAGGAAGATAGGCATGTAGCTTGATTACCTACATGTCTATCGCCAGGGAAATG 27426 
ck/CH/LDL/091022 TT T- AG 27474 


Fig. 2. ( Continued ). 


organizations to those of other IBVs. In a phylogenetic 
analysis using the full-length genomic sequences, both ck/ 
CH/LZJ/111113 and ck/CH/LDL/091022 were clustered into 
the same group (Fig. IB). 

3.3. Recombinant analyses of ck/CH/LZJ/111113 genome 

Because no full-length genomic sequence of 4/91-type 
was available, we pairwise compared the nsp 13, S, M, N 
and 3' UTR of the 4/91 strain with those of ck/CH/LZJ/ 
111113 and ck/CH/LDL/091022. As shown in Fig. 2, the nsp 
13 of the ck/CH/LZJ/111113 strain bears a much greater 
nucleotide sequence identity to that of ck/CH/LDL/091022 
(96.3%) than to 4/91, though the ck/CH/LZJ/111113 1 -1453 
in the SI gene showed an identical sequence to that of 4/ 
91. It worth noting that in agreement to the nsp 13, strain 
ck/CH/LZJ/111113 1453-3492 in spike, M and 1-1015 in 
the N genes shared 99.3, 99.0 and 94.9% identities to those 
of ck/CH/LDL/091022, respectively, in contrast to the 84.8, 


91.3 and 86.9% to 4/91. To our surprise, ck/CH/LZJ/111113 
from 1032 to the end of the N gene and 3' UTR again have 
identical sequences to those of 4/91 and were different 
from those of ck/CH/LDL/091022 (91.3 and 75.2%, respec¬ 
tively). In addition, ck/CH/LZJ/1 111 13 and ck/CH/LDL/ 
091022 had a 96% nucleotide identity; however, this figure 
was up to 97.2% when the 4/91-like sequences in the 5' 
ends of SI and N were excluded. These data strongly 
suggested that ck/CH/LZJ/111113 arose from a homolo¬ 
gous RNA recombinant event from multiple template 
switches (at least three switches) between 4/91- and ck/ 
CH/LDL/091022-like viruses. Here, we investigated the 
presence of two switch sites at the 3' ends of the SI and N 
genes, respectively (Fig. 2), implying that template 
switches occurred within the SI and N genes. Due to the 
unavailability of the sequences between nsp 13 and the 
spike genes in the 4/91 virus, we could not identify the 
exact switch site before the spike gene; however, by 
comparing the sequence with that of ck/CH/LDL/091022, 
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ck/CH/LZJ/i 11113 nsp 13 gene 


ck/CH/LZJ/111113 

100 

ck/CH/LDL/091022 



ck/CH/LDL/091022 

96.3 

100 

4/91 

4/91 

88.1 

89.6 

100 


ck/CH/LZJ/111113 


spike gene 

ck/CH/LZJ/111113 

100 

ck/CH/LDL/091022 



ck/CH/LDL/091022 

90.8 

100 

4/91 

4/91 

91.2 

82.5 

100 


ck/CH/LZJ/111113 

S fragment before switch site 

ck/CH/LZJ/111113 

100 

ck/CH/LDL/091022 



ck/CH/LDL/091022 

78.7 

100 

4/91 

4/91 

100 

78.7 

100 


ck/CH/LZJ/111113 


S fragment after switch site 

ck/CH/LZJ/111113 

100 

ck/CH/LDL/091022 



ck/CH/LDL/091022 

99.3 

100 

4/91 

4/91 

84.8 

85.1 

100 


ck/CH/LZJ/111113 


M gene 

ck/CH/LZJ/111113 

100 


ck/CH/LDL/091022 



ck/CH/LDL/091022 

99.0 

100 

4/91 

4/91 

91.3 

91.2 

100 


ck/CH/LZJ/111113 


N gene 

ck/CH/LZJ/111113 

100 


ck/CH/LDL/091022 



ck/CH/LDL/091022 

94.5 

100 

4/91 

4/91 

88.9 

87.3 

100 


ck/CH/LZJ/111113 


N fragment before switch site 

ck/CH/LZJ/111113 

100 


ck/CH/LDL/091022 



ck/CH/LDL/091022 

94.9 

100 

4/91 

4/91 

86.9 

86.5 

100 


ck/CH/LZJ/111113 


N fragment after switch site 

ck/CH/LZJ/111113 

100 


ck/CH/LDL/091022 



ck/CH/LDL/091022 

91.3 

100 

4/91 

4/91 

100 

91.3 

100 


ck/CH/LZJ/111113 


3 untranslation region 

ck/CH/LZJ/111113 

100 


ck/CH/LDL/091022 



ck/CH/LDL/091022 

75.2 

100 

4/91 

4/91 

100 

75.2 

100 


Fig. 2. ( Continued ). 
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we could predict that this site may tend to lie immediately 
upstream of the S glycoprotein gene because ck/CH/LZJ/ 
111113 and ck/CH/LDL/091022 shared highly identical 
sequences in this region. 

3.4. Virulence studies 

Clinical signs were observed in all ck/CH/LDL/091022- 
infected chicks from days 3-14 post-inoculation and in all 
ck/CH/LZJ/111113-infected chicks from days 4 to 18 post¬ 
inoculation. The clinical signs included listlessness, hud¬ 
dling, ruffled feathers and dark, shrunken combs. The 
clinical signs were more severe in ck/CH/LZJ/1 111 13- 
infected chicks than those infected with ck/CH/LDL/ 
091022. Two of the 10 chicks died: one on day 4 and a 
second on day 7, in the ck/CH/LDL/091022-infected group 
during the experiment; however, no chicks died in ck/CH/ 
LZJ/111113-infected group. Gross lesions in the organs of 
the dead chicks were confined mainly to the kidneys. No 
overt disease was observed in chicks of the control group. 
In addition, chickens inoculated with ck/CH/LZJ/111113 
and ck/CH/LDL/091022 were negative for IBV antibodies 
until 8 days post-inoculation. At each time point from 8 
days post-inoculation onwards, the ck/CH/LZJ/111113- 
inoculated chickens had comparable antibody titers with 
those of ck/CH/LDL/091022. 

3.5. Analysis of IBV-RNA in tissues of inoculated chickens by 
real-time PCR 

As illustrated in Fig. 3, all tissues investigated had 
detectable viral RNA on days 4 and 7 post-inoculation, 
respectively. Overall, the replication efficiency of ck/CH/ 
LDL/091022 in the chicken tissues was higher than that of 
ck/CH/LZJ/111113. However, both strains showed a some¬ 
what high affinity to the gastrointestinal tract, as reflected 


by the high viral genomic copies detected by real-time RT- 
PCR. Furthermore, ck/CH/LDL/091022 had a strong affinity 
to the kidney, compared to that of ck/CH/LZJ/111113. No 
IBV-derived RNA was detected in the tissues of the control 
chickens. 

4. Discussion 

In 2011, we conducted routine surveillance of IBV in 
China and isolated and identified an IBV strain from an 
H120-vaccinated chicken, designed as ck/CH/LZJ/111113. 
Serological data suggested that ck/CH/LZJ/111113 pos¬ 
sessed distinct characteristics setting it apart from the 
predominant LX4-type strain in China (Liu and Kong, 2004) 
and the mass-type vaccine strain. Phylogenetic analysis of 
the ck/CH/LZJ/111113 strain showed that it was mostly 
related to pathogenic UK 4/91 and the recently isolated 
Indian virus strain India/NMK/72/IVRI/10, which strongly 
suggested that this isolate was of the 4/91 type. However, 
the clearly divergent 3' end of the SI gene, compared to 
that of 4/91, implicated that a recombinant event may have 
occurred during the origination and evolution of the virus. 
Further proof of multiple recombinant events was found 
by extensive analysis of the full-length genomic sequences. 

Coronaviruses are known to efficiently recombine, 
especially within the spike gene, to modulate pathogen¬ 
esis and even the host range. This has been well 
documented in IBV, which was clearly shown to be 
responsible for the emergence of several new serotypes or 
variants (Jackwood et al., 2010), and contributed to the 
outbreaks of diseases in chickens. The recombination 
events from which the ck/CH/LZJ/111113 virus resulted 
can be explained by two models. In the most simple 
scenario, the recombination may have involved only two 
parental viral strains with RNA replication initiating in a 
ck/CH/LDL/091022-like template of either negative or 
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Fig. 3. Viral genome RNA copies in tissues from chickens inoculated with ck/CH/LZJ/111113 and ck/CH/LDL/091022, respectively. Viral RNA copy number 
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positive polarity and then switching of the polymerase- 
nascent cRNA complex to a 4/91-like virus template, 
followed by double switches at the 3' ends of the SI and N 
genes, respectively. The first switch may have occurred in 
the vicinity of the ORF 1/S boundary, presumptively 
upstream of the spike glycoprotein gene, but the exact site 
for this is currently unknown because of the sequence of 4/ 
91 at this position is unavailable. In this scenario, the ck/ 
CH/LDL/091022- and 4/91-like viruses must be both 
circulating in the animal breeding facilities for long 
periods, in which chickens were infected simultaneously 
and recombination events occurred. Alternatively, a more 
complicated scheme can be envisaged. In this scheme, a 
ck/CH/LDL/091022-4/91 hybrid arose from a single switch 
at either of the three or from two switches at two sites and 
spread among the chicken flocks and in turn engaged in 
additional recombination events with another 4/91-like 
virus. IBV circulates extensively in chicken flocks, either 
broiler or layer hens, and even in the vaccinated flocks in 
many countries, in some cases, with more than one type of 
strain. Conceivably, this and the long persistence of IBV in 
chickens could raise the odds of occurrences of double or 
triple infections. 

It is believed that the conditions for recombination 
amongst the IBV strains in the field are as follows: an 
extremely large number of chickens, most maintained at 
high density; ease of spread of the virus; and serotype co¬ 
circulation, including proof of co-infection with more than 
one serotype in a given flock. This is true in the case of ck/ 
CH/LZJ/111113. China has an extremely large number of 
chickens and most of the chickens in different districts in 
China are maintained at a high density. LX4-type IBV 
isolates are found to be the predominant IBV type 
circulating in both laying hens and broilers in China. 
Though only small numbers of 4/91-like viruses have been 
isolated in chicken flocks in China (Han et al., 2011), 
however, 4/91 live vaccine are commonly used during the 
breeding period, the vaccines must certainly be implicated 
as a source for 4/91-like sequences. The LX4- or 4/91-like 
sources of genetic material were probably a strain that 
naturally infects these affected flocks, as found in ck/CH/ 
LZJ/111113 in this study. The natural recombination event 
demonstrated here could suggest that the continuous use 
of live vaccines may actually contribute to natural 
recombination and IBV-associated disease. 

Overall, the body of evidence in favor of a recombinant 
origin for the SI, partial N genes and 3' UTR of ck/CH/LDL/ 
091022, together with its distinct antigenic and tropism 
profiles, argue in favor of setting it apart from other known 
IBVs as novel strains. However, a recombination event 
itself is unlikely to be the only reason directly responsible 
for these distinctions between ck/CH/LZJ/111113 and ck/ 
CH/LDL/091022 because only a few mutations were found 
in other parts of the genomes. Because the S glycoprotein 
gene is involved in host cell attachment and contains viral 
neutralizing epitopes, recombination of the S glycoprotein 
gene can result in the emergence of new strains or 
serotypes of the virus as well as new viruses with distinct 
characteristics such as pathogenicity and tissue tropism. 
This may likely to be the actual reason for the change in 


serotype and tissue tropism for the ck/CH/LZJ/111113 
strain compared to ck/CH/LDL/091022. 

In the CoV genome, as most RNA viruses, the 5' and 3' 
UTRs usually harbor important structural elements that 
are involved in replication and/or translation (Masters, 
2006). In IBV, the 3' UTR binds to the N protein which is 
essential for synthesis of negative-strain viral RNA. 
Perhaps the acquisition of the 3' end of the N gene and 
the 3' UTR from 4/91 by ck/CH/LZJ/111113 can alter the 
viral replication efficiency, as observed in this study. 
Hence, this alternation may in turn affect pathogenicity. 
However, further investigations using reverse genetic 
system should provide further insight into this issue and 
increase our understanding of IBV pathogenesis. 
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